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We find the propagator and calculate the tree level scattering amplitude between two 
covariantly conserved sources in an Anti-de Sitter background for the most general D- 
dimensional quadratic, four-derivative, gravity with a Pauli-Fierz mass. We also calculate 
. the Newtonian potential for various limits of the theory in flat space. We show how the 

' recently introduced three dimensional New Massive Gravity is uniquely singled out among 

. higher derivative models as a (tree level) unitary model and that its Newtonian limit is 

' equivalent to that of the usual massive gravity in flat space. 
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^ ; I. INTRODUCTION 

In gravity, there seems to be an insurmountable difficulty in reconciling renormalizability with 
^IIh . unitarity in generic dimensions. By adding higher derivative terms, aB? + PRfiu^ to the four 

I dimensional Einstein-Hilbert action, one gains renormalizability yet looses unitarity due to a non- 

CN ' decoupling ghost introduced by the /? term, without which one does not have a perturbatively 

. renormahzable theory The efi'ect of non- unitarity in the Newtonian hmit shows itself as a 

Q I repulsive component to gravitational force between static sources. Because of this repulsive com- 

' ponent, the theory has a better UV behavior. This is what usually happens in field theory: ghosts 

. are introduced to make the theory better behaved, at least during the process of renormalization, 

I yet, in a unitary theory they decouple at the end. Of course, bartering unitarity with renormaliz- 

Q>^ ' ability, as in the case of four dimensions, cannot be accepted. 

■ In three dimensions the situation seems to be better as was recently shown in 0! S]: with 

the choice 8a + 3(3 = and a reversed sign Einstein-Hilbert term, one obtains a perturbatively 
renormahzable, 'unitary' theory in fiat space [1, [E)Sj0]- But, it is not clear at all if this particular 
^ ■ ratio between a and /? will survive renormalization at a given loop level, even at one-loop. The 

, unitarity beyond tree level has to be checked. What is also interesting is that the linearized 

version of the theory has in its spectrum a massive graviton with helicities ±2. This fact sheds 
light, albeit only in three dimensions, to an old problem of finding a non-linear extension to the 
Pauli-Fierz mass term. In fact, a formal equivalence of the Einstein-Hilbert-Pauli-Fierz gravity 
and the linearized version of the "New massive Gravity (NMG)" was shown in 11]. [Note that, 
in three dimensions besides this parity-preserving spin-2 theory, we have the old parity violating 
Topological^ Massive Gravity (TMG) [3,0].] Even though there does not seem to be a problem 
with this formal equivalence, one still needs to be careful about its physical meaning. It is clear that 
even the linearized version of the NMG theory is background diffeomorphism invariant, but the 
Pauli-Fierz theory is only invariant under the Killing symmetries of the spacetime (in particular, 
the 2-|-l dimensional Minkowski space). Therefore, a better understanding of the symmetries is 
needed. A quite interesting approach was put forward in where, in the linearized version of 
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the NMG theory without the Einstein-Hilbert term, Weyl invariance of the action was shown. 
Therefore, at least in the hnearized level, introduction of the Einstein-Hilbert term breaks this 
invariance and introduces a mass to the graviton. From this point of view, higher derivative terms 
provide the kinetic energy and the Einstein-Hilbert term provides the mass in this model, which 
also explains the bizarre sign change of the Einstein-Hilbert action. In retrospect, this is to be 
expected, pure Einstein's theory is non-dynamical and gives no propagation in three dimensions: 
at the linearized level, it is basically like the mass term in a scalar field theory which can only 
play role in the dynamics once a kinetic energy is introduced. Here the kinetic energy comes from 
the higher derivative terms. So both in TMG and NMG, Einstein-Hilbert term gives rise to the 
mass in the linearized theories, by breaking Weyl invariance ( not the expected diffeomorphism 
invariance). This point of view could be important for constructing massive gravity theories in 
other dimensions. 

After [3], several works appeared that were devoted to different aspects of this theory. Especially, 
in AdS background that we shall be interested in, some classical solutions and conserved charges 
were presented in 0, H, E E Q. In this paper, we consider the most general quadratic 



model, augmented with a Pauli-Fierz mass term, in a -D dimensional (anti)-de Sitter background. 
We study the propagator structure and the tree level scattering between two covariantly conserved 
sources in the theory obtained by linearizing the action 

/ = y d^x - ^ + + pjf^^ + ^ (^2^^^ _ 4^2^ ^ I 
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+ d'^x^gl-— {hi, - h') + £inatter > (1) 



4k 

where Aq is the bare cosmological constant and k is related to the D— dimensional Newton's 
constant and the D — 2 dimensional solid angle by k = 2Q.£,^2Gd- Including the number of 
dimensions, this 7-parameter theory is the most general, four derivative, quadratic model with 
various potentially interesting limits and discontinuities. [In the absence of the source terms and 



at the linearized level, one can reduce the number of parameters in the action |l5l . Il6l |. but here 
for the sake of generality we shall work with ([1]).] At this point, we assume nothing about the 
signs of the parameters in the action, moreover we will also allow them to vanish. Constraints will 
come from the requirement of the tree- level unitarity and the non-existence of ghosts and tachyons. 
For some specific dimensions, certain terms will not contribute to the equations, for example, the 
7 term (Gauss-Bonnet combination) is a total divergence in = 4 and vanishes identically for 
= 3, therefore, in three dimensions, the Riemann tensor carries no more information than the 
Ricci tensor. Also in Z) = 3, one can add the Chern-Simons term ^i{TdT + f F^) g, Q to extend 



our model, but here we will stick to ([T|), since TMG is special to three dimensions. [See 171. Il^ 
for this case, without the higher curvature terms.] For D = 2, the theory reduces to an model 
with a Pauli-Fierz term. Here, we will consider D > 3. Apart from the Pauli-Fierz mass term, the 
theory has general covariance. 

The spin-2 model defined by the linearization of ([1]) is highly non-trivial. Needless to say, 
various limits have been studied in the literature. Yet, there still appears interesting new models 
in certain limits of the above action. NMG, for D = 3, in the case of a flat background, and for 
(M^ = 0, 7 = 0) and 8a -|- 3/3 = 0, being one such example. One of our tasks in this paper is 
to explore in detail the full 7-parameter theory with particular care on the various discontinuities 
that appear in changing the order of limits when some of the parameters approach zero. One very 
well-known discontinuity is the so called van Dam-Veltman-Zakharov (vDVZ) discontinuity: the 
fact that Einstein-Hilbert gravity in flat space is isolated from massive gravity in a discontinuous 
manner, that is 0, at the tree level does not yield the correct General Relativity (GR) 

results. But once a cosmological constant is introduced and M^/A — > limit is taken, GR result 
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is recovered 19, [2ol . 21, 22, [2^. [Another resolution of the discontinuity may fohow even in flat 
space if the Schwarzschild radius of the scattering objects is taking as a second mass scale in the 
theory 23|].] But, these are ah at tree level, once quantum corrections are taken into account, 
discontinuity reappears 2j]. Another related problem is the Boulware-Deser instability: at the 
non- linear level, a ghost arises in massive gravity [25]. 

The lay out of the paper is as follows: In section II, we write down the linearized equations 
about an AdS background and discuss certain special limits, such as the partially massless case. 
In section III, we find the tree-level scattering amplitude and compute the Newtonian potential. 



II. LINEARIZED EQUATIONS 



The field equations that follow from ([T]) are 

^ (^R^u - ]^g^^R + ^og^,^ + 2aR [r^^ - ^g^,uR^ + (2a + 13) {g^^U - V^V^) R 



+27 



RR. 



'^RaavpR^'' + R, 



pcfpr 



R, 



crpT 



2R 



ml. + R^) 



+I3U (r^. 



-g^,R ) + 2/3 



R 



M2 



gpuh) 



' flU ■ 



(2) 



In the absence of the source, the 'vacuum' ( background ^^jy) is a non-singular solution to the field 
equations everywhere. This is the maximally symmetric (anti)-de Sitter Space with the Riemann, 
Ricci tensors and scalar curvature given respectively as 

2A 5 2A _ - 2DA 

[gpugpa - gfiugvp) , R^u — — ^^g^u, R — -f^ — -■ (3) 



R 



ppi/a 



Our discussion and notations follow 2^, 27 1. All the contractions will be made with g^i^ that has the 
signature (—,+,...,+). Our conventions are [V^,Vt/]Vx = Rp.ux'^y^rr, Rpu = R'^ p.\v With these, 
one can find A in terms of a,f3,'y,K and Aq (neglecting M^), where A is the effective cosmological 
constant. From ^ and ([3]) we have. 



A 



2k 



D 



(Da + p) D 
+ 7- 



A^ = 0. 



- ^. (4) 
D -2 ' D -1 ^ ^ 

Generically, there are two solutions one of which vanishes in the absence of Aq, the other one 
nonvanishing even in this limit. But, lets first note the exceptions. For D = 4, A = Aq. Also, for 

D = 3 if 3a + /? = we have A = Aq. Otherwise, in three dimensions Ai^2 = ^^^^^^g'^^j^.'^j"^^^" - with 
a constraint 1 > 8k (3a + (3) Aq. There is one more exceptional point: 7 = 0, Da + /3 = theory 
also has A = Aq. [This is an interesting model whose action is given by the square of the traceless 
Ricci tensor. All the asymptotically AdS solutions ( not just the globally AdS vacuum) have zero 
energy 271] .] Finally, for the generic case we have 

1 



A 



4k/ 



1 ± Vl + 8k/Ao 



(5) 



where / 

Linearization of 



{aD + (3) ^^-2 + 7 |^_'^||^_2) 1 for which the bound becomes as SkAq/ > 



{D-2r 

2|) around the background metric, g^^i, 



-1. 



Tf,u (h) = agj:^ + (2a + /3) g^^^O - V^V^ + 



gp.u + h^i,, after using gives [271] 
2A 



D 



^gpu 



R' 



L 
pu 



2A 



D 



jgpi^R 



+ 



M2 

2k 



{hp.u - gp.uh) , 



(6) 
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where we have defined 



1 AKD 4A ^ 4Am-3)(D-4) 

— I a H d H — 

K D-2 D-l^ {D-l){D-2) ^ 



(7) 



Here T^^ (h) contains all the higher order terms as well as the source t^^. Q^v^ is the linearization 
of the Einstein tensor (with the cosmological constant) 



R 



1 



2A 



L 



where the linearized Ricci tensor and the scalar curvature, R 



R 



{g'"'R^^ur, read 



2A 



D 



(8) 



■h. (9) 



We also need the trace of (161) which reads 



{Aa{D - 1) + D/3) □ - (D - 2) ( - + 4/A 



R' 



{D - l)h = 2T. 



(10) 



It is clear that something special happens for 4a{D—l)+D(3 (that is the NMG point in D = 3 ) but, 
before, we discuss this in detail by computing the tree level scattering between two sources, let us 
consider various limits of the theory at the linearized level, without the sources. Unitarity regions 
cannot be captured this way, but we can see the parameter ranges that rule out the tachyons. At 
this point, the discussion depends on whether vanishes or not. Let us consider these cases 
separately. 

/ case 



Divergence and the double divergence of ([6]) give 



V - = 0, 



Dh = 0, 



(11) 



which lead to R'" = —^^h. The discussion again bifurcates according A = or not. First consider 
the flat space case for which we have R^ = and (jlOp and the first equation of (llip force the field 
to be traceless and transverse. The field equation for the remaining {D + 1){D — 2)/2 independent 
components reduces to 



K K 



0, 



that describes two massive excitations with masses 



1 



m 



± 



2kP 2\kP\ 



\/1 + 4/3M2k, 



(12) 



(13) 



which are non-tachyonic if the parameters are properly picked. But, as we will see, this model is 
nonunitary. 

Now assume that A 7^ 0, then the trace equation gives 



1 \ 

iia{D - 1) + D(3)n - (D - 2) ( - + iAf\ + ^{D -1){D -2) 



h = 0. 



(14) 



which says that generically /i is a dynamical scalar field, unless the coefficient of the D'Lambertian 
vanishes. Suppose, we pick up that special case, then, there are still two options: either h = 
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or we have the partially massless point, that arises only in curved backgrounds for which a 
higher derivative gauge invariance appears [28] and the field has one less degree of freedom (DOF) 
compared to the massive one. The mass should be tuned as 

2 2Ak fl Al3{D-A) _ ^^_{D -3){D -A) 



which is allowed to be negative in AdS as long as it satisfies the Breitenlohner-Freedman type 
bound. Apart from four dimensions, higher derivative terms play role on the partially massless 
theory. 

= case 

The theory is now invariant under background diffeomorphisms 5^hfj_y = V^^^u + VuC/j., since 
S^Qj^iy = and 6^R^ = 0. Therefore, divergence and the double divergence do not give any 
constraint on h^,y. In this case, for T = 0, (jlOp gives dynamics to unless, the coefficient of 
the box term vanishes. In that special case, generically = 0, but it is clear that cosmological 
constant introduces another possibility : If ^ + 4A/ = 0, then R^ need not vanish. But this point 
does not seem acceptable since a gauge invariant object is left undetermined by the field equations. 



III. TREE-LEVEL AMPLITUDE 

From now on, we consider the full theory ([6]) and find the tree level scattering amplitude between 
two covariantly conserved sources. First we need to express hfj^iy in terms of Tfj_^. But, since not all 
components of the field are independent, we decompose it in such a way that the physical parts 
will be determined by the source. The usual choice is to define 

V = ^lu + V(^K) + v^v^</. + Qf^^^P, (16) 

where /i^J is the transverse and traceless part. Symmetrization (with a 1/2 factor ) is implied 
in the vector part which is divergence free, (j) and ip are scalar functions. Taking the trace, 
divergence and double divergence of Eq. (fT6l) one obtains 

2A - - 

h = D(l) + Dtl;, nh = &(!)+ 50 + 5-^, (17) 

[L> 2) 

where we used V'^V^/i^jy = Oh, which is not a gauge condition but imposed on us as a result of 
the nonzero mass term. Then, hitting the first equation of ()17p with a □, one can eliminate 
with the help of the second equation as 

p^^ (»-l)(»-2) p^_ (18) 



which then yields 



/.h '^'T'" "^"!^- (19) 



From (jlOp . it follows that "0 is determined by the trace of the energy momentum tensor 
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where c = ^^^_^2°' + "U^- ^^'^ transverse traceless part of the field in terms of the source, 

(2) 

Lichnerowicz operator, A]^ acting on spin-2 symmetric tensors is quite useful: 

V = -QV - '^Rf^puahf'" + 2RP (21) 
Some properties of this operator that we need were collected in [i^] 

A?V(^K) = V(^a(;V,), a(;V^ = (-□ + A) V^, V^Af V = Ai'V'^ V> 

Af = 9^.a(°V, aJ^V = -a<P, V^A^JV^ = A^PV^^y^. (22) 

Using these we have 

f Ail/ ~ 2 ^ (L) - 2) '^^ ' 

which then leads to the desired equation 

= 2 {(/?□ + a)(A(^) - ^) + TJJ. (24) 

The transverse traceless part of the energy momentum tensor can be found as ( after using the 
fact that it is covariantly conserved ) 



D-1 D-iy " {D-l){D-2) 



2AD ^ 



Finally using (j20l 1241 [25]) . we can write the tree level scattering amplitude between two conserved 
sources 

A = \jd^x ^T'^, (x) h^^^ i^) = \J d^'x ^ (t;,/^^^^- + TV) . (26) 
For the sake of notational simplicity, lets suppress the integral for now: 

iA = 2r;,|(/3n + a)(Af --^) + ^| 

4A „, 4A . M^l ^ 2AL> ^ 



r { {pa + a){n + ——) } ^5+ . — -} T 



(D - 2)(L> - 1)2 {''^ D-2' kJ t {D-2){D-l 

2 r (i + 4A/ - - - 1)1" (□ + ,„ .V TXT,) 



{D-2){D-1) [Hi ' 2kA' 'j { {D-2){D-l 

This is our main result from which we will consider various limits. For nonzero cosmological 
constant, this is quite a nontrivial integral. But we can figure out the particle spectrum of the 
theory by looking at the pole structure of the amplitude. Generically there are 4 poles which read 
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as 



°^ - ~{D-l)^D-2y ^^^^ 
□2,3 = ^<;-(^^ + 7#4)±l/(^^ + 77=^) -/^(ttF^-— )^ (29) 




□4 = _U-i + 4A/-— . (30) 



Given these poles, finding the residues is easy but, in the most general form the expressions are 
rather cumbersome. The existence of the cosmological constant and the Pauli-Fierz mass term 
changes the picture drastically: Depending on the choice of the parameters, there could be tachyon 
and ghost-free models. Here, let us restrict our theory and discuss some interesting limits and 
compute the Newtonian potential between static sources for some of these limits. 

Looking at \27\ . it is clear that — > and A ^ limits do not commute. In fact first taking 
the A ^ limit, one encounters the vDVZ discontinuity. 

AA = -2T'{[5d^ + -d^-—\ + -l-T'[pd^ + -d^ -—\ \, (31) 

f^" y K kJ D-l [ K K ] ^ ' 

whose spectrum has two massive excitations with masses given by p3p . To see its structure more 
explicitly, one can rewrite it as 

4.A = (32) 

f3im'i - ml) \ - ml - ml J (D-l) yd'^ - ml d'^-mlj 

Unless /3 = 0, we have a massive ghost. The Newtonian potential energy (U) between Tqq = 
mi6{x — xi), T^^ = m26{x — X2) in three and four dimensions can be obtained as 

U = -pP^\K^(m^r) - K,{m^T)\ D = 3, 

2p(ml — m_) 47r 

U = —^f^^^-Lle-""-^ - e-'^+n D = 4. (33) 

3p{ml - ml) Airr^ ^ ^ ' 

where r = \x\ — X2\- As /? ^ 0, the potential energies become 

U = -—mim2Ko{Mr) D = 3, (34) 



AGm,m2 Mr ^ ^ (35) 

3 r ^ ^ 

The latter equation shows the famous discontinuity of massive gravity in fiat space. [Note that 
we have used k = IGvrG, in four dimensions.] Let us stress that (j34p is the Newtonian limit of 
massive gravity in three dimensions. It gives an attractive force as long k is positive. Unlike the 
four dimensional case, M ^ limit does not exist since, as x ^ 0, Kq{x) — > — ln(x/2) +je, which 
incidentally gives the expected 1/r force for small separation between the sources. So massive 
gravity gives the correct Newtonian limit in three dimensions where pure Einstein theory does not 
give any interaction. 
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Let us now consider the case which wih lead us to the NMG model found in [5]. Take first 
= then A ^ 0. 



-1 



2 T'\cd^--d^} T (36) 



{D-l){D-2) { K 
Generically there are three poles : 

,1 = 0. a| = -i .81 = 1^ (37) 

For nonzero /?, the requirement of unitarity singles out the D = 2> and 8a + 3/3 = theory in 
a highly nontrivial way. Let us explain how: not to have a tachyon, we should choose Kf3 < 0. 
Looking at the residue of this pole, we see that we should have k < for unitarity. Computing 
the residue of the massless pole, we see that, with negative k, D = 3 must be chosen so that one 
does not have a massless ghost. Finally, the third pole is non-tachyonic if c < 0, but then, the 
residue of this pole requires c > for unitarity. This is only possible if c = 8a + 3/3 = 0, which is 
the NMG theory. Newtonian limit of this model also reveals its rather unique structure. Let us 
for the moment compute the potential for generic a and /3. 

U = ■^mim2 {Ko{mgr) - Ko{mor)) D = 3, (38) 

ovr 

where rrig = — ^ and ttiq = ^(gQ,+3/3) • Clearly, mo is a massive ghost that gives a repulsive 
component. But, for NMG it decouples and one is left with an attractive force, since k < 0. This 
result also confirms that, at this level, NGM has the same Newtonian limit as the usual massive 
gravity p4p . if the Pauli-Fierz mass term is chosen as M = nig. Beyond three dimensions, in flat 
space, massive ghost does not decouple unless /? = 0. As an example, let us look at D = 4: 

U = (l - y-^^ + ) , (39) 

where = 2K(3a+)3) • '^^^ niiddle, repulsive term signals the ghost problem 

IV. CONCLUSION 

We have studied the most general quadratic gravity with a Pauli-Fierz mass in D dimensional 
(anti)-de Sitter space. From the tree level scattering amplitude that we found, one can study 
various limits. In flat space, we computed the Newtonian limits for various models including the 
New Massive Gravity that was recently introduced. Non-unitarity of the NMG theory comes in a 
highly non-trivial way and does not extend beyond three dimensions, in flat space. The cosmological 
constant changes the picture drastically, one needs to further study in detail the unitary regions. 
Especially in the NMG theory, as we mentioned in the introduction, unitarity beyond the tree level 
has to be checked, it is not clear at all if the condition 8a + 3(3 = will survive renormalization, 
nor it is clear that Boulware-Deser instability in the full non-linear theory is avoided. We intend 
to address these problems in a future work. 
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